A Magnetic Model for the Ordered Double Perovskites 
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We construct an effective spin model from the coupled spin-fermion problem appropriate to double 
perovskites of the form A2BB'06- The magnetic model that emerges is reminiscent of double 
exchange and we illustrate this 'reduction' in detail for the case of perfect B-B' structural order, 
i.e, no antisite disorder. We estimate the effective exchange between the magnetic B ions in terms 
of the electronic parameters, study the 'classical' magnetic model using Monte Carlo techniques, 
and compare this approach to a full numerical solution of the spin-fermion problem. The agreement 
is reasonable, and promises a quick estimate of magnetic properties when coupled with ab initio 
electronic structure. The scheme generalises to the presence of antisite disorder. 



I. INTRODUCTION 

Double perovskite materials, of the form A2BB'06, 
have been of interest in recent yearsi^ on account of 
their magnetic, electronic and structural propeties. They 
promise large magnetoresistance^-, potentially useful 
for switching applications. The half- metallic character of 
some of the members also make them attractive candi- 
dates for spintronic devices. 

One of the species, B say, is typically magnetic, a tran- 
sition metal like Fe, Co, Ni, or Cr, while the B' species is 
generally non-magnetic, Mo, W, etc. The most studied 
member of this series is S^FeMoOe: it is a half-metallic 
ferromagnet (FM) at low temperature, and has a high 
T c ~ 410K. Si^FeWOe, on the other hand, is an anti- 
ferromagnetic (AFM) insulator! These limits illustrate 
the wide range of physical properties in the double per- 
ovskites (DP). While the 'endpoints' above are relatively 
easy to understand (ignoring disorder) there are several 
effects where current understanding is limited. 

(i) Antisite disorder: Well annealed double perovskites 
tend to have an alternate arrangement of B and B' ions, 
but defects called 'antisite' regions^ appear when two B 
or two B' atoms occur as neighbours. These regions typ- 
ically have an AFM arrangement of the B spins and are 
insulating. Their presence reduces the overall magneti- 
zation. The electronic and magnetic properties in DP's 
are intimately related to the structural order. 

(ii) Phase competition: Exploration of the series^ 
Sr 2 FeMoi_ a; W a ;06 reveals a FM to AFM transition and 
an associated metal-insulator transition (MIT) with in- 
creasing x. In the regime of FM-AFM phase competition 
the compounds show large magnetoresistance (MR). 

(hi) Magnetic B' sites: Recently, compounds where the 
B' site also has an intrinsic magnetic moment have been 
investigated^, and interesting compensation effects have 
been observed. In particular, there are enigmatic com- 
pounds like S^CrOsOe which are insulating (semimetal- 
lic?), but at the same time ferromagnetic, with a very 
high T c 10 . In addition, there are spin-orbit effects^ in 
some DP's complicating the magnetic state. 

Apporaching issues (i)-(iii) above directly in a hnite 
temperature real space formulation is formidable. It re- 



quires tools that can predict magnetic properties of a 
double perovskite based on electronic parameters and 
the structural disorder. This paper is a step towards 
that goal where we provide a semi-analytic scheme for 
accessing the magnetic ground state and T c scales of 
a structurally ordered DP starting with a tight-binding 
spin-fermion model. While our primary focus is the FM 
regime, we also highlight issues of phase competition and 
antiferromagnetism which are bound to be important 
when doping effects are explored. 

The paper is organised as follows. The next section 
describes the double perovskite model, following which 
we summarise earlier work on this problem to place our 
work in context. We then outline the different meth- 
ods used in this study. The section after describes our 
results, primarily within a variational scheme and an ef- 
fective exchange calculation, with Monte Carlo results 
for benchmark. We then conclude, pointing out how our 
scheme can be extended to the antisite disordered case. 



II. THE DOUBLE PEROVSKITE MODEL 

The double perovskite structure of A2BB'06 can be 
viewed as repetition of the perovskite units ABO3 and 
AB'03. In the ideal ordered DP the B and B' octahedra 
alternate in each direction. In this paper we consider only 
the B ion to be magnetic. The superexchange coupling 
between the B magnetic moments is small in the ordered 
DP's. The important physical ingredients in the problem 
are: (i) a large S core spin at the B site, (ii) strong cou- 
pling on the B site between the core spin and the itinerant 
electron, strongly prefering one spin polarisation of the 
itinerant electron, and (iii) delocalisation of the itinerant 
electron on the B-O-B' network. 

The Hamiltonian for the structurally (B-B') ordered 
double perovskites is given by: 



i£B i£B' i 

- t Y fl m ^ + J J2 Si ■ fhvtpfif} (1) 

<ij>a i£A 

The /'s refer to the magnetic B sites and the m to the non 
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magnetic B', and the B-B' hopping tBB 1 = t is the prin- 
cipal hopping in the structurally ordered DP's. We will 
discuss the impact of further neighbour hoppings later in 
the text. We have retained only one orbital on the B and 
B' site, our formulation readily generalises to a multiple 
orbital situation. The Si are 'classical' (large S) core 
spins at the B site, coupled to the itinerant B electrons 
through a coupling J> i. This implies that the conduc- 
tion electron state at a B site is slaved to the orientation 
of the corresponding B spin. The difference between the 
ionic levels, A = e B — e B ,, defines the 'bare' 'charge trans- 
fer' energy. At a later stage we will define the parameter 
A = {(.b — JS/2) — 6b> as the 'true' charge transfer en- 
ergy, n/ is the B electron occupation number, while n m 
is the B' electron occupation number. We will assume 
J/t — > oo (keeping A finite). The parameter space of 
the problem is defined by the electron filling, n, the ratio 
A/t and the temperature T/t. We have ignored Hubbard 
repulsion, B-B antifcrromagnetic superexchange, and, to 
start with, direct hopping between B'-B', or B-B. 



III. EARLIER WORK 

Early work on the DP's was motivated by results on 
Sr2FeMoC>6, where electronic structure calculations indi- 
cate that Fe is in a 3c? 5 configuration (a half-filled state) 
while Mo is in a 4d 1 configuration. Following Hund's rule, 
Fe is therefore in a high spin S = 5/2 state. Surprisingly, 
the normally nonmagnetic Mo picks up a moment of 1/2 
in the opposite direction, and reduces the moment per 
unit cell to ~ 4/xg. An explanation for the induced mo- 
ment on the non magnetic B' species was provided^ by 
Sarma et al, in terms of a 'level repulsion' between the 
Fe and Mo levels. Such a scenario implies a substantial 
degree of hybridization between the Fe and Mo orbitals, 
and assumes that the itinerant Mo electron hops through 
the Fe sublattice. 

Using this idea, a double-exchange (DE) like 2- 
sublattice Kondo lattice model was proposed for 
the DP'si^, and solved within dynamical mean field 
theory by Chattopadhyay and Millis^, assuming a 
'ferrimagnetic— state. They obtained a n — T phase 
diagram for different values of cb> — e_B and J and ob- 
served that the T c goes to zero at large filling, indicating 
the presence of some competing non ferromagnetic state. 

A similar result for T c (n) was obtained by Carvajal et 
ofiS, using another two sublattice model, and Ising spins. 
Here the hopping of an electron with spin a from a B' site 
to a neighbouring B site is i if a is antiparallel to the local 
spin Hi on that site, while it is zero if they are parallel. 
The authors considered only ferrimagnetic arrangements. 

Alonso et a££ considered a variant of Millis' model 
with the coupling J — * oo, but with a larger number of 
ordering possibilities. They also took into account pos- 
sible antisite defects, including a B-B hopping and su- 
perexchange which are only active when two B atoms 
become nearest neighbours. They considered four possi- 



ble phases: (1) paramagnetic, (2) ferrimagnetic, (3) an 
AFM phase, where the B spins in neighbouring (1,1,1) 
planes are antiparallel, and (4) another ferrimagnetic 
phase where the B spins are aligned ferromagnetically if 
the B are in the correct positions, and antiferromagneti- 
cally if the B ions occupy B' sites due to antisite defects. 
Among other results they found that even in the B-B' 
ordered case (where superexchange is not operative) the 
AFM phase is preferred to the FM at high band filling. 

All these studies, except the paper by Alonso et al., 
concentrate on the ferromagnetic^ phase. They observe 
the decrease of T c at large filling but do not explore com- 
peting phases. Secondly, while the DMFT approaches 
provide a semianalytic treatment of the T c scales, in spe- 
cific parts (in this case ferromagnetic) of the phase dia- 
gram, an estimate of the effective exchange between the B 
moments is not available. The ab initio approaches have 
attempted such an estimate by force fitting a 'Heisenberg 
model'. Unfortunately, the magnetic states that emerge 
from the DP model, and the effective exchange that sta- 
bilises these phases arise from subtle electron delocalisa- 
tion physics not captured by such methods. We also do 
not know of any work that allows an economnical and 
systematic exploration of the parameter space, n, A, J, 
of the DP model. The present paper aims to overcome 
these shortcomings. 



IV. METHODS 

The first estimate of magnetic interactions in any ma- 
terial is provided by ab initio calculations. This is 
typically done by calculating the difference in ground 
state energy of the compound in spin polarized and 
spin-unpolarizcd configurations; or in different magnetic 
ground states corresponding to different values of the spin 
density wave vector—. Such a calculation involves all 
the relevant orbitals and their hybridization and provides 
a rough material specific estimate. However, for com- 
plex antiferromagnetic ground states one has to guess 
such configuration beforehand, or take a cue from ex- 
periments. There is no a-priori prescription for finding 
them. 

Model Hamiltonian based calculations, on the other 
hand, have the obvious limitation that model param- 
eters have to be inferred from elsewhere, typically ab 
initio studie s 18 ' 19 . The advantage, however, lies in the 
simplicity of the resulting model, and our ability to 
create a qualitative understanding using the tools of 
statistical mechanics. The Hamiltonian appropriate to 
double perovskites can be studied using the following 
tools: (i) a combination of exact diagonalization and 
Monte Carlo (ED-MC), (ii) variational calculation (VC) 
based on some family of periodic spin configurations, and 
(iii) mapping to an effective classical spin model. 

The ED-MC approach has the advantage of accessing 
the magnetic structure without bias. However, due to 
large computational cost, it is severely size- limited, limit- 
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ing the class of magnetic structures which can be probed. 
A 'travelling cluster' (TCA) variant 2 ^ of ED-MC allows 
use of somewhat larger system size. Variational calcula- 
tions assuming a periodic spin background can be used 
for very large system size (since there is no bulk diago- 
nalisation needed) but are restricted by the choice of the 
variational family. While we will use both (i) and (ii) 
above, our principal tool will be (iii), where we map on 
the spin-fermion problem to an effective spin only model, 
with exchange calculated from the fermions 2 ^. We de- 
scribe (i)-(iii) in more detail below. 

A. Monte Carlo 

One can solve the DP model on a finite lattice by direct 
numerical methods, allowing for an 'exact' benchmark 
for approximate solutions. ED-MC is such a technique. 
Here, the coupled spin-fermion problem is solved by up- 
dating the classical spins using a Monte Carlo, diagonal- 
izing the fermion system at each step of the MC to infer 
the energy cost of the move. The method is numerically 
expensive and can only be used on small system sizes, 
~ 8 x 8. Substantially bigger sizes, ~ 24 x 24, can be 
accessed using the TCA. 

For the MC implementation the Hamiltonian of Eq [1] 
has to be cast into form appropriate for J — > oo. This 
is done by performing a rotation to the local Si axis at 
each B site, and retaining only the electron state oriented 
antiparallel to Si at that site. This gives the following 
Hamiltonian, with 'spinless' B conduction electrons and 
B' electrons having both spin states. 

H = t £ {(«n(|)/tm iT - e^cos(^)fj mil ) + h.c.} 
<ij> 

i ia 

There is no longer any 'infinite' coupling in the model, 
and the number of degrees of freedom has been reduced 
to one per B site (and 2 per B'), so the Hilbcrt space is a 
little smaller. vn,j\ and hop to different conduction 
electron projections at the neighbouring B site(s) so the 
effective hopping picks up a 9i , fa dependent modulation. 
We will use this form of the DP model for the Monte 
Carlo. 



B. Variational ground state 

A more analytical method used before in the double 
exchange context is to write down a family of spin con- 
figurations {S} Q , denoted S a for simplicity, and calculate 
the electronic energy in that background. Since the S a 
are usually periodic this is effectively a 'band structure' 
calculation. For a specified chemical potential one can 
calculate the electronic energy £(fi, S a ). The configura- 
tion S rn i n (fi) that minimises £ is the variational ground 



state. Needless to say, the 'minimum' is only as good as 
the starting set, and in general non periodic S a cannot 
be handled. Nevertheless, used in combination with MC 
results it can be a valuable tool. 

From the MC we will discover that in the structurally 
ordered case the DP model has simple periodic ground 
states, with windows of phase separation in between. 
This will allow us to use the variational scheme, with 
only a few configurations, to map out the T = phase 
diagram accurately. 

C. Effective exchange 

The complications with spin-fermion MC, and the lim- 
itations of VC could be avoided if one had an explicit 
spin-spin interaction model deduced from the starting 
DP model. Formally such a scheme can be written down, 
and some progress made through approximation. Let 
us illustrate this 'self consistent renormalisation' (SCR) 
principle 26 in the simpler context of double exchange be- 
fore moving to the double perovskites. 

1. Illustrative case: double exchange model 
Consider the following model: 

H = *« C L<> - J S * • AcfiapCiP (3) 

Let us try to construct an approximate classical spin 
model in the limit J — * oo. The classical model is de- 
fined by the equivalence: 

j DSie-^'W = J VSiTre-? 11 (4) 

where the trace is over the fermion degrees of freedom. 
The trace, in general, is impossible to compute analyt- 
ically since it involves the spectrum of fermions moving 
in an arbitrary spin background {S}. Nevertheless, some 
headway can be made once the Hamiltonian is written in 
a more suggestive rotated and projected basis as 2 ^: 

H = Y,fi3UA^ ii ih + M (5) 

ij 



where — y — ^ — 1 , <£>i j is a phase factor depending 
on Si and S,-, and the 7 are 'spinless' fermion operators. 
This suggests the approximation: 

H eff {S} « -E^V 1± f^ 

Da = -t i j((e i *»'rh+h.e)) (6) 

The angular brackets indicate first a quantum average 
(for fixed {S}) and then thermal average over e~^ He ff^ . 
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Another way to obtain the same result, which gener- 
alises to the DP problem, is to write the action for H in 
a spin background {S}: 

.4{S} = p V [i^Sij - iy/ye**] 4 n l jn (7) 



and the internal energy U{S} = dln ^ s ^ = _(cL4) ; 

U{S} = E&^^E^n) 

ij n 

= E/yttf^frh^ + M (8) 

which is simply the quantum average of the spin-fermion 
Hamiltonian for a fixed {S}. As before we can convert 
this to an approximate spin Hamiltonian by thermally 
averaging the quantity within the round brackets. 

The effective exchange depends on T, but in the 'clean' 
problem it does not depend on the 'bond' ij. Since the 
DE model always has a ferromagnetic ground state, the 
low T exchange can be calculated from the fermionic av- 
erage in the fully polarised state, and is simply: 

D oc ^£fe(n fc ) = y^^e k n F (e k ) 

k k 

It was observed^ that even at finite temperature the self 
consistent average in the ferromagnetic phase remains 
close to the T = value till very near T c . The T = 
kinetic energy therefore provides a reasonable estimate 
of effective ferromagnetic exchange, and so the T c . The 
overall scale factor between the T c and the exchange can 
be determined from a Monte Carlo calculation. 



where tk = 2t^2 i=l coskia. The poles of this Greens 
function give the band dispersion at J = 0: 



± A±^ 

h k - 7T 



4e 



/'• 



(12) 



In the limit A>4, i.e., the limit of weak charge transfer, 
there are two bands centred roughly on and A. For 
A = 0, there are two bands ±|ejfe| symmetrically placed 
about 0. 

While the first term in the action involving this bare 
Greens function conserves spin and momentum, the sec- 
ond term is local in real space and typically involves spin- 
flip. To proceed, let us Fourier transform Gjj Q (k, iS) and 
write the action in real space. t\ generates 'hoppings' 
(in the full B-B' lattice) connecting sites that can either 
be next nearest neighbours (2N), next-to next nearest 
neighbours (3N), or the same site. In real space the ac- 
tion assumes the form: 

A{S} = ^(EfLGj/ott-^^fs™ 

J^Si-finaVapfinp) (13) 



n %j a 



Now, an unitary transformation is performed in spin 
space so that the second term in the action becomes di- 
agonal: 7i njU = J2 a A ^afina- The action becomes: 

-4{S} = J2 9%,-rt ltn Gj} {f i -i f j ,iu> n )~( jvn 



JS 



-ltl n liln)) (14) 



2. Effective exchange in the double perovskites 

Unlike the DE model we cannot write an effective spin 
only Hamiltonian for the double perovskites purely by 
inspection since the electron motion also involves the B' 
sites. We use the action formulation instead. Integrating 
out the B' electrons we get an action entirely in terms of 
the B degrees of freedom: 

,4{S} = /?£(£/L^7/VM^)/fe« CT 



n ko 



J E S i • fin^vfinv) (9) 



where Gffo(k,iuj n ) is the J = Greens function involv- 
ing B sites only (the n represent Matsubara frequencies): 



//0 " " t ~' iu n ~ (es' - M) 

If we choose es = 0, eb> = A < 0, this becomes: 



Gf f0 = iuj n + jU - 



iuj n + [I - A 



(10) 



(11) 



where o y = V A 1 A 3 

At large J one projects out the ji U states, retaining 
only the terms involving the index I. Thereafter, we drop 
this index, redefine the B level as eb — > es — 4r, and 
obtain an effective spinless fermion model similar to the 
case of double exchange. 

MS] = P E 9ijJin G Jfo( f i ~ ? 3 ' iuJ n)l 3 n 
ij 

foi. i 



ico„ + [i — A 



)lvnljn 



where flj,- = 



as before, ha is the Fourier 



transform of t\ and connects sites on the B sublattice. It 
involves a NN term (x + y in the full B-B' lattice) and a 
third neighbour term (2x etc in the B-B' lattice). 

It is important we appreciate the various terms in the 
expression for ^4{S} above. The 'kernel' Gjf (ri—fj,ioj n ) 
is specified by the J = bandstructure of the B-B' prob- 
lem, explicit information about the spin variables is en- 
coded in g t j, and the fermions are defined in the back- 
ground {S} . 



FIG. 1: Colour online: The three magnetic phases in the 2D model. Left: ferromagnetic (FM), center: antiferromagnet 
(AFMf), right: antiferromagnet (AFM2). These occur with increasing electron density. The moments are on the B sites, we 
have not shown the induced moments on the B' sites. 



Wc define iv n = i(2n + 1)tt, so 



The internal energy can be calculated as J7{S} = 
~^§JT = ~( 9 "a/? S * )- Simplifying the resulting expres- 
sion and using the same principle as in DE we can write 
an explicit (but approximate) model purely in terms of 
core spins: 



Heff{S} 

B(iw n ) 



£a. 



1 + Si.Sj 



B(iu n ){e i9 » (7 in 7jn) + h.c) 

n 

{2iu n + 2fi - A) 
(iu n + fi- A) 2 



(15) 



The effective exchange can be determined at any 
temperature by the SCR method. The couplings take 
two values, D\ for NN B-B exchange, and D2 for sec- 
ond neighbour B-B exchange. The low T exchange can 
be estimated by evaluating the fermionic average in the 
perfectly spin-ordered state at T = (after checking 
that the ground state generated by the exchange is self- 
consistently ferromagnetic). The evaluation of the Mat- 
subara sums, etc, is discussed in Appendix A. 



V. RESULTS 

A. The magnetic ground state 

Both ED-MC done on 8 x 8, and TCA done on 16 x 16 
exhibit the presence of three phases: namely ferromag- 
netic (FM, first panel in Fig.l), a 'line like' antiferromag- 
nctic phase (AFM1, middle panel) and the more con- 
ventional antiferromagnet (AFM2) in the last panel. If 
we define the ordering wave-vector on B sublattice us- 
ing axes along the diagonals, the FM phase has order at 
Q = {0,0}, AFM1 has order at Q = {0,tt}, and AFM2 



has order at Q = {n,ir}. We could of course define the 
wave- vectors on the full B-B' lattice and use the usual 
x and y axes, but for the structurally ordered case the 
earlier convention is simpler. 

Fig [5] shows the n — T phase diagram for three values 
of A. With increasing n the phases occur in the sequence 
FM, AFM1, AFM2, AFM1 and FM again. The sequence 
as well as the rough filling windows are similar for all 
three A values (the VC, which is free of size effects, will 
demonstrate this more clearly) . Fig IHb) shows that for 
intermediate level difference, A = 4, the T c for the ferro- 
magnetic phases actually increase a little bit. The T c of 
the antiferromagnetic phases, however, decrease. More- 
over, the {0, 7r} phase is unobservable on the high filling 
side, while its Neel temperature, Tjy, is quite small even 
on the low filling side. Eventually, for large enough A, 
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FIG. 2: Colour online: n — T phase diagram based on TCA. 
(a) A = 0, (b) A = -4, and (c) A = -10. The solid (blue) 
lines are ferro, dashed (red) lines are AFM1 and the dash-dot 
(black) lines are AFM2. The system size is 16 x 16. 
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the T c of even the ferromagnetic phases decrease, as seen 
in Fig[^c). The AFM1 phase is unobservable on even 
the low filling side, possibly due to very small T^- 

While AFM phases driven by B-B superexchange have 
been studied in the DP's, AFM phases driven by electron 
delocalisation have not seen much discussion. Their oc- 
curence, however, is not surprising. If we were to 'test 
out' the feasibility of various magnetic ground states we 
could restrict ourself to a few simple collinear phases 
to start with. The FM, AFM1, etc, are such exam- 
ples. Let us index them by some index a. As de- 
scribed before, which of these occur at a chemical poten- 
tial \x can be simply checked by calculating the energy 
£<*(m) = f-oodcN a (e)e, where N a (e) is the electronic 
density of states in the spin background a. The phase 
appropriate to a particular [i would be the one with low- 
est energy. Even without a calculation it is obvious that 
the FM state will have the largest bandwidth, and would 
be preferred at low n. The AFM phases have narrower 
bands, but larger density of states (since the overall DOS 
is normalised), and with growing /i they become viable. 
In what follows we quantify this carefully. 

One can obtain analytic expressions for the dispersions 
in the {0, tt} and {tt, 7r} phases, which are given below. 
In the {0, 7r} phase, in our J — > oo limit, the structure 
decomposes into electronically decoupled ferromagnetic 
zigzag chains aligned antiferromagnetically with respect 
to each other, see Fig.l middle panel. Their dispersion 
is ID-like, given by: 



A ± J A 2 + 16 + l6cos(k x - k y ) 
£ k = - ^ (16) 



In the limit A — ► 0, their ID like nature is clearly visible: 

(17) 



e k = 2v / 2c OS (^— ^) 



The {n, n} phase, on the other hand, decouples into 
two planar lattices where the B spins are arranged ferro- 
magnetically, while these lattices are themselves aligned 




FIG. 3: Colour online: Electronic density of states for the 
three variational states, with Q = {0, 0}, {0, n} and {tt, tv} 
at A = 
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FIG. 4: Colour online: Electronic energy (top) and filling 
n (bottom) versus chemical potential for the three phases, 
A = (left), A = -2 (right) 



antiferromagnetically with respect to each other. The 
dispersion is given by: 



£k = 



A ± y/A 2 + l6t 2 (cos 2 k x + cos 2 k y ) 



(18) 



cos 2 k y when A = 0. It 



which reduces to ±2ty // cos 2 k x 
is interesting to note that the bandwidths of both {0, tt} 
and {it, tt} phases are identical, although the detailed 
DOS are different. The DOS for the three phases for 
A = are shown in Fig [31 We can understand the oc- 
curence of the various phases by integrating the DOS 
and comparing the energies at a fixed /i. The results are 
shown in Figg^a) for A = and Figg^b) for A = -2. 

As expected, at low filling the energy of the FM phase 
is the lowest, while for intermediate filling that of the 
AFM1 phase is lower than the FM phase. At still higher 
fillings, the energy of the AFM2 phase is the lowest. This 
is repeated symmetrically on the other side of [i = 
for A = 0. The density discontinuity corresponding to 
each transition can also be found from the corresponding 
fj, — n curves. For finite A, the AFM1 phase becomes 
narrower, especially on the high filling side. These simple 
variational results are corroborated by the phase diagram 
obtained from the TCA calculation. 

The n — A phase diagram at T = is shown in Fig [5] 
There are windows of phase separation (PS) where ho- 
mogeneous electronic/magnetic states are not allowed. 
These regions correspond to the jumps in the n— [i curve. 
The AFM1 phase becomes unstable on the high filling 
side for large A, which manifests itself through a merg- 
ing of the phase boundaries. 
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FIG. 5: Colour online: n — A phase diagram at T — from the 
variational calculation. The hashed regions indicate windows 
of phase separation between the adjoining phases. 



B. Spin model and effective exchange 

Our effective spin model is: 

ij 

Using the exchange calculated from Eq [15] using 
a fully ferromagnetic reference state, one can plot the 
nearest neighbour exchange Dnn and the next nearest 
neighbour Data? at as a function of filling n. The results 
are shown in the top panel of FiglBJa) for A = 0. One 
finds that both the exchanges change sign as a function 
of filling. At low filling and very high filling, both are 
negative, indicating an overall ferromagnetic coupling. 
However, for intermediate values of filling, both the ex- 
changes become positive, giving an effective antiferro- 
magnetic coupling. In between, there is a small region 
where one of them is positive and the other negative. 

Since the calculation was started using a purely ferro- 
magnetic spin background, such changes in sign of the 
calculated exchange indicate an instability of the ferro- 
magnetic phase at these fillings. Where the 'exchange' 
Dnn+Dnnn > the ground state will no longer be FM, 
the result for is not self-consistent, and the quanti- 
tative values not trustworthy. We will confine ourself to 
the window where the ground state is self consistently 
ferromagnetic. 

The exchange for the antiferromagnetic states, and the 
Neel temperature, should be calculated in appropriate 
spin backgrounds, i.e., {0,7r} and {tt, tt}. However, the 
{0, 7r} state, in the J — ► oo case, consists of disconnected 
chain-like structures. While the intra-chain arrangement 
is ferromagnetic, the inter-chain arrangement is antiferro- 
magnetic. Since there is no hopping connectivity between 
the chains, the inter-chain exchange calculated in such 
a spin background would emerge to be zero. Similarly, 
for a {tt, tt} spin background, there are two sublattices 
such that the intra-sublattice arrangement is ferromag- 
netic, while the inter-sublattice one is anti-ferromagnetic. 
Again, since these sublattices are disconnected, the inter- 
sublattice hopping is zero. In these anisotropic states the 




FIG. 6: Colour online: Aj for NN and NNN, A = (top) 
and n — T phase diagram obtained from a Monte Carlo on 
the spin model using this exchange (bottom). 



effective exchange (and stiffness) vanishes along certain 
directions at T = 0. To calculate the effective exchange 
that controls the T c in the AFM1 and AFM2 phases 
we need to necessarily solve the finite temperature self- 
consistency problem. This is an interesting problem, but 
computationally demanding, and is left for future work. 



C. n — T phase diagram from the Dij 

Our previous experience with the double exchange 
model suggests^ that a reasonable estimate of T c is pro- 
vided by the exchange calculated in the fully FM T = 
state. On this assumption, one can study the effective 
spin model with classical Monte Carlo and calculate finite 
temperature properties including T c . The n — T phase 
diagram obtained this way is shown in the bottom panel 
of FigH All the three phases: FM, AFM1, and AFM2 
occur in approximately the correct filling windows. The 
T c scales for the ferromagnetic phases, which are the only 
ones consistent with the assumed spin background, turn 
out to be reasonably correct, as we will see in a compar- 
ison with the full TCA result. We ignore the T c for the 
AF phase since the AF exchange is not self-consistent. 



D. Properties in the ferromagnetic regime 

Since much of the interest in the double perovskites 
arises from ferromagnetism, we focus on this regime in 
what follows. In our n — T phase diagram, this FM phase 
at low filling occurs upto n « 0.5 — 0.7. From the varia- 
tional calculation, which is essentially in the 'bulk limit', 
the FM window is upto ~ 0.3 — 0.4. Considering the 
degeneracy of the three t 2g orbitals, translates to about 
0.9 — 1.2 electrons per unit cell. S^FeMoOg, which has 
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FIG. 7: Colour online: Ferromagnetic part of the phase dia- 
gram for low 'filling' 17 , x — 3 — n, compared between TCA 
and SCR, for (a) A = 0, (b) A = — 10 both calculated for size 
16 x 16 




FIG. 8: Colour online: Effective exchange in the ferro phase 
\Dnnn + Dnnnn\ for different A, obtained from SCR. The 
exchange calculation is on a k grid 1000 x 1000, and the 
fillingii is x = 3 — n. 



one electron per unit cell, falls within this regime. How- 
ever, many other materials like Si^FeReOg are known, 
which have 2 or more electrons per unit cell, but are still 
ferromagnetic with a high T c . This discrepancy between 
theory and experiment was noticed by many authors be- 
fore us: Chattopadhyay and Millie 13 , L.Brey et. a/. 20 , 
E. Carvajal et.al^ and J.L.Alonso et. al*£. They at- 
tributed it to the presence of competing antiferromag- 
netic channels, an effect which we also find. 

Our results on T c (n) is similar to that obtained by 
others, i.e, a reduction as n — > as the kinetic energy and 
ferromagnetic exchange weakens, and a drop also at large 
n due to the presence of competing AF phases. The T c (n) 
obtained from the SCR scheme is compared to the result 
of full spin-fermion Monte Carlo using TCA, Fig[7J They 
seem to match quite well, except that the SCR results 
calculated on a T = state overestimate the T c slightly. 



While the actual T c -s can only be calculated using Monte 
Carlo on small systems, eg. for 16 x 16 in Fig the 
average exchange can calculated with a very large k grid 
(lOOOx 1000 k points). This T = exchange, based on an 
assumed FM state has the behaviour shownii in Fig [8] 
While this result has a clear correspondence with the 
T c (n,A) results obtained from DMFT by Millis et. al, 
and the calculations of Carvajal et.ai, it overestimates 
the window of FM, and misses the first order FM to AFM 
transition. The actual T c (x) will have a discontinuity 
with increasing x, instead of decreasing smoothly to zero. 

E. The AFM phases 

The AFM phases AFM1 and AFM2 occupy a large 
part of the n — T phase diagram. The presence of such 
collinear antiferromagnetic phases have been observed 
earlier by several authors, notably Alonso et.al. 16 . These 
phases, at least within the J — ► oo model considered here, 
have a a 'lower connectedness' than the ferromagnetic 
phase. The AFM1 phase consists of double staircase-like 
structures attached back to back, while the AFM2 phase 
consists of decoupled Cu-0 like lattices for each B' spin 
channel. The DOS corresponding to these are given in 
Fig O The DOS for the AFM1 phase resembles that of 
a ID tightbinding lattice, while the AFM2 DOS is more 
2D-like. It is interesting to note that there is a dispersion- 
less level for both the AFM phases, which gives the jump 
in the fi — n curve. While the effective exchange calcula- 
tion starting from a fully polarized background already 
produced the three phases, but a truly self-consistent cal- 
culation for the AFM phases would have to start assum- 
ing these spin backgrounds. Such a calculation is non- 
trivial, as discussed before. 

VI. MAGNETISATION AT THE B' SITES 

The induced magnetism on the B' site was explained 
within a local 'level repulsion' picture by Sarma et a&^, 
but a lattice-oriented approach is lacking. Some headway 
can be made by exactly integrating out the B degrees of 
freedom rather than the B' from the J — > oo model given 
in Eq[2l but the result is in the form of an action, rather 
than an effective Hamiltonian. Within second order per- 
turbation theory, however, it appears that an extra onsite 

term of the form YljSafi S j+S ' mj a u aP m iP occurs 

for the B' sites (z is the number of nearest neighbours) 
giving an 'exchange splitting' at the B' site, with the 
spins of the surrounding B sites serving as a magnetic 
field. The corresponding hopping terms are given by: 

W / t , ^FM a t - 

<ij>,<? <ij> 

t 2 

+ 2A ^ S ^i)/ 2 ■ m L a ^m jf3 (19) 

<<ij>> 
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FIG. 9: Colour online: Temperature dependence of B' site 
magnetization for different filling, from ED-MC on a 8 x 8 
system at A = 0. The top panel plots the magnetisation 
normalized by the total B' occupancy, while the bottom panel 
shows the unnormalized magnetisation. 




0.025 0.05 0.075 

T 



FIG. 10: Colour online: Comparison of the electronic mag- 
netisation at the B' site with the core spin magnetization at 
the B site. The result is for n = 0.4 using ED-MC on a 8 x 8 
system. Both the B and B' results are normalised to highlight 
their similar temperature dependence. 

Within a mean-field treatment of the B core spins Si w 
(Si) = zM, the effective B'-B' 'hamiltonian' can be writ- 
ten as: 

t 2 

- [(1 + M)m^m n + (1 - M)ml imjl ] 

<ij> 
t 2 

+ ^ [(! + M K ! T m iT + (1 - ^KW] 

<<ij» 

2t 2 

+-^ M ^2( m it m n - m n m ii) ( 20 ) 

i 



Obviously, at T = 0, M = 1, and only one spin species 
hops. 

The effective spin polarization at the B' site, which is 
purely electronic, contributes to the total magnetization. 
Within the TCA approach using the Hamiltonian [21 this 
can simply be estimated by calculating the normalized 

magnetization — 57 57 — . In Fig [5Ja) , the magneti- 

Th .j. ~|" Th 1 

zation of the B' has been plotted against the temperature 
for different fillings correponding to an ED-MC simula- 
tion on a 8X8 system for A = 0, normalized by the net B' 
filling. It is observed that the T dependence is very sim- 
ilar to that the B core spin case. In Fig[9](b), the bare B' 
magnetization is shown without normalization: it shows 
that the saturation magnetization increases with filling, 
as expected. In FiglTUl a comparison of the M vs T com- 
ing from the B' electron and the B core spins is provided. 



VII. THE EFFECT OF B'-B' HOPPING 

Inclusion of B'-B' hopping t! would result in the same 
expression for the exchange calculated in the ferromag- 
netic state as before (see Appendix A), except that the 
A everywhere would get replaced by A + At' cosk x cosk y , 
while ek = 2t(cosk x + cosk v ) for a square lattice. It is 
obvious that if t = 0, the exchange would vanish irrespec- 
tive of t', showing that hopping across the magnetic site 
is crucial, as expected. For parameter values reasonable 
in double perovskites, t' w 0.1 — 0.3t, and < A < 3, 
there is almost no change in the n — A phase diagram, 
although the T c values decrease marginally when the t' is 
turned on. For larger values of A, the AFM1 phase be- 
comes unstable, and the ferromagnetic window extends 
a bit, upto the AFM2 phase, although the T c s, of course, 
are proportionately low. 



VIII. DISCUSSION 

We discuss a few issues below to connect our results to 
available data on the double perovskites, and also high- 
light a few effects that we have neglected. 

1. Material parameters: Ab initio calculations 
sugg est±2ii2 that t ~ 0.3 - 0.5eV, while t' is typically 3-5 
times smaller. The direct hopping between B sites is even 
smaller, ~ 0.05eV. Estimates for the bare 'charge trans- 
fer gap' A (in SFMO) vary between 1.4eV— to about 
2eV— . Hence, the parameter window we explored seems 
reasonable. Our ferromagnetic T c are typically O.li at a 
filling appropriate to SFMO, so the absolute magnitude 
of the T c -s would be about 360 — 600A', roughly the range 
seen in the for double perovskites. 

2. B-B hopping: While we have not considered the 
effect of B-B hopping, the smallest energy scale in the 
problem, explicitly in this paper, it is possible to under- 
stand qualitatively the effect of including this hopping. 
If only B-B' hopping is considered, then there are two 



10 



singular features in the density of states, at es and cb', 
or alternatively, at and A. Inclusion of the B'-B' hop- 
ping resulted, at the zeroth level, in the smoothening 
of the feature at cb 1 - Similarly, inclusion of the B-B 
hopping will basically smoothen out the feature at €b- 
However, the B'-B' hopping had a much more dramatic 
consequence in terms of providing an alternate pathway 
for derealization of the B' electrons irrespective of spin, 
especially at large A. Secondly, it resulted in connecting 
up of the AFM1 staircases, getting rid of their ID char- 
acter, and making this phase unstable compared to the 
FM and AFM2 phases. Inclusion of a small B-B hopping 
in addition, is not, on the other hand, expected to have 
any more dramatic consequences. We can readily include 
this in our formalism. 

3. Three dimensions: The entire analysis in this pa- 
per was in two dimensions. Apart from simplicity, ease 
of visualization, and computational tractability, there is 
a definite argument in terms of the symmetry of the ti g 
orbitals as long as nearest neighbour interactions are con- 
sidered—^, which says that one can consider three inde- 
pendent 2D Hamiltonians. Other author s 21 ' 23 have also 
used 2D Hamiltonians. Having said that, the phases dis- 
cussed here generalizes easily to three dimensions. The 
AFM1 phase in the absence of B'-B' hopping becomes 2D 
rather than ID, consisting of ferromagnetic [111] planes 
arranged antiferromagnetically. Such an arrangement 
has been observed by authors like Alonso et. al^, and 
even in ab initio calculations — . 

4. Filling control: The n — T phase diagram that 
we provide is for a definite set of parameters A, t, t' 
etc. While going across the series of DP compounds 
Sr 2 FeMo06, to Sr 2 FeRe06, it is not just the filling but 
also all these parameters which are changing. A more 
controlled way of varying the filling alone would proba- 
bly be to dope the compounds at the A site, namely pre- 
pare the series S^-zLa^FeMoOg. While some work has 
been done in this regard 2 ^, more extensive work, probing 
higher doping values is necessary to ascertain whether 
such antiferromagnetic phases are indeed observed. 

5. SCR for antisite disordered case: While the clean 
problem has been studied in detail in this paper, antisite 
disorder is expected to make it more interesting. Para- 
doxical effects like increase in T c and widening of the 
FM region has been suggested 16 . The scheme for self- 
consistent renormalization that we have proposed can be 
generalized even to the case of antisite disorder (see Ap- 
pendix) . However, the scheme in that case becomes more 
numerical, and the analytical handle available here would 
be lost even at T — 0. The formalism is presented in the 
Appendix and we are currently studying the problem. 



IX. CONCLUSION 



The 'exchange scale' in this model is related to the elec- 
tronic kinetic energy. The ferromagnetic T c estimated 
from this exchange compares well with results from the 
full spin-fermion Monte Carlo. The change in sign of 
the exchange with increasing carrier density captures the 
phase competition in the electronic model and the transi- 
tion in the magnetic ground. Our scheme, extended to in- 
clude multiple bands and spin-orbit coupling would allow 
a controlled and economical approach to the finite tem- 
perature physics of a wide variety of double perovskites. 
Another fruitful line of exploration is to include anti- 
site disorder. We have highlighted the scheme in the 
appendix and hope to present results in the near future. 

We acknowledge discussions with D. D. Sarma, Bri- 
jesh Kumar and Rajarshi Tiwari, and use of the Beowulf 
cluster at HRI. 



X. APPENDIX A: EXCHANGE CALCULATION 
IN THE DOUBLE PEROVSKITES 

The effective exchange can be evaluated in the per- 
fectly spin-ordered state at T=0. This can be obtained 
by using the known form of the Green's function (7 i „7 J n) 
at T=0, namely it is the Gjfo(fi — fj, iu> n ) obtained be- 
fore, made dimensionless by dividing by (3. Hence, in the 
spin ordered case, the exchange becomes: 



(3 ^ (iw. 



{2iuj n - A) 



(2iu n - A)e 2 



-y- 

f3 ^ (lUJ n - A)[lLO n {lUJ r , 



(21) 



Expanding in partial fractions, this can be written as: 

Ek+ E k+ A 



U B 



-y 



P 



iuj n - EkA 



iuj n - E k - 



- A 



Performing the Matsubara sums, this gives the final re- 
sult 

U B = y [E k+ n F {E k+ ) + E k _n F {E k _)} - An* (A) 
k 

The last term is an additional contribution obtained 
from the missing B energy. While this gives the full 
internal energy in the spin polarised case, the bond- 
resolved exchange Dij is given by 



Di 



^VM^-r;) [E k+nF (E k+ ) + E k _n F (E k _)} 

k 

— An F (A) (24) 



We have suggested a scheme for extracting a simple 
magnetic model for double perovskites starting with a 
tight binding parametrisation of the electronic structure. 



We need to do the k-sum, which can only be done 
numerically for square (2D) or cubic (3D) lattices. In- 
stead, if one uses a Bethe lattice of infinite coordination, 
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then the bare density of states for nearest neighbour hop- 
ping is semicircular, and analytic treatment is possible, 
at least for the DOS. In our case, there are two bands, 
with dispersions 



Ek± 



A ± y/A 2 + 4e 2 



(25) 



The density of states for these bands is: 



J duipo{oj)5 ( E 



A ± VA 2 + 4w 2 



(26) 



where the bare DOS po(u>) = ^2 k S(lo — ek) is approxi- 



mated by the semicircular DOS, po(u>) ~ 
Then, the full DOS given by Eq |2"6"1 becomes: 

1 !' u , —, -,(,-. A ± %/A 2 + 4w 2 ' 



P±{E) 



ttD 2 



J duy / D 2 -cj 2 5 (e 



(27) 

This integral can be evaluated, to give the following an- 
alytic result for the DOS: 



P±{E) 



±(2E - A) y/D 2 - E' 2 + EA 



ttD 2 



VE 2 - EA 



x 0{D 2 - E 2 - EA)G(E 2 - EA) 



(28) 



Obviously, this diverges at E = 0, and E = A. The limits 
of the DOS, and hence the integral, are obtained from 
the equations given by the two theta function conditions: 
= Effik = 0, A If we take A < 0, as 

in our case, then the lower band lies between A +^A 2 +4 i> 
(left edge) and A (right edge), while the upper band lies 
between (left edge) and (right edge). 

When D — > 0, i.e., t — ► 0, the theta function conditions 
are only satisfied together for E 2 — EA = 0, i.e., at E = 
or E = A. Thus we recover the bare levels as (5-function 
peaks in the DOS, as expected. 

The exchange at T=0 is given in terms of this DOS as: 

U B = J E[p + (E)+p_(E)]n F (E)dE-An F (A) (29) 

This gives Ub as a function of fi. One can also obtain the 
total number of electrons ./V from the DOS as a function 
of /x: 

N = J dE[p+(E) + p_{E)]n F (E) (30) 

From Eq |2"91 and Eq[3Dl eliminating p, one can get U vs 
N. 



Using the substitution u> — E 2 — EA, the expression 
for the exchange can be rewritten in a more convenient 
form at T=0 as: 



Ub(p) 



A ± VA 2 + Auj 2 



x© [p 



A ± VA 2 + Aoj 2 



VD 2 



AQ(p - A) 



(31) 



It is to be noticed that as t — > 0, i.e., D — > 0, the ex- 
change goes to zero, as it should. This can be seen in two 
ways. Firstly, as t — > 0, e^. — > 0. Hence, the band disper- 
sions Ek± given by Eq [25] tends to A= ^ A ^ . This means 
that E k+ — > and E k _ — > A. Hence, putting in Eq |2"91 
the first term involving E k+ is 0, while the second term 
involving E k - cancels the third term involving An F (A). 

The other way to observe this is to use Eq[2U Here, 
when D — > 0, then the theta function condition is only 
satisfied for uo 2 — D 2 , which means that the only con- 
tribution to the integral comes from lo = 0. Indeed, the 
bare semicircular DOS -^jyiVE) 2 — uj 2 being a normalized 
object, tends to a delta function 5(E) as D — > 0. Hence, 
the term involving + sign gives 0, while that involving — 
sign gives A6(^- A), which cancels with the third term. 



XI. APPENDIX B: EXCHANGE 
CALCULATION WITH ANTISITE DISORDER 

Firstly, the Hamiltonian has to be written in such a 
way that all the B and B' degrees of freedom are sepa- 
rated out in distinct subspaces of the Hamiltonian. 



H = 



Hpp Hpm 
Hmf Hmm 



where Hpp represents the terms in the subspace of the 
B degrees of freedom, while Hmm represents the terms 
in the B' subspace. Hmf and Hpm connects the two 
subspaces. 

The B Green's function satisfies the matrix equation 

G FF (iuj n ) = iuj n I — Hpp — HpM^rJ- — Hmm) -1 Hmf 

(32) 

Written out term by term, 



G FF (iu> n ) — iui n S i j—Hpp ij —'^^HpM ik (iunI—HMM)ki 1 HMF i . 

hi 

(33) 

Hence the action 

in ij 

~f3 2 ^ Hp M . k {lV n I - [3H M m)Zi HMF zj linljn9ij 
kl 
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Now, let V be the diagonalizing matrix of Hmm and 
its eigenvalues are A^ fM . Then, 



(j > j 



in ijn 

HFM ik Vk p V pl 1 Hn 



/dA 

kip 



-^E 



l MF u 



kip 



9ijlinlj n 



£ 

ij 



H 



FF, 



H F AU k V kp V t 1 H M f 1j 



j / j 

kip 



H F M ik V kp \ P M M V pl 1 H M F h 



AIM 



(w n - px p MM y 



9ij < linljn > 



It is to be noted that in the ordered case, the hamilto- 
nian matrix becomes block diagonal in k-space, and the 
elements of the off-diagonal block Hmf and Hfm & re 
simply 6^, while those of the digonalizing matrices V are 
e ik.(n-rj) ^ wn ii e t ne eigenvalues \ V MM are simply fj/o, 
i.e., A. Hence, the quantity HFM ik Vk p V p ^ Hmf^ in the 
numerator simply goes over to hij, as defined in Eq 1151 
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